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It  is  currently  common  practice  in  theoretical  ocean  acoustics  to  derive  parameter  resolution  bounds 
for  a  monochromatic  measurement  of  the  temporally  fluctuating  field  received  by  a  hydrophone 
array.  However,  a  monochromatic  measurement  corresponds  to  a  single  random  sample.  In  applied 
ocean  acoustics,  single  samples  are  seldom  if  ever  used  for  parameter  estimation  because  the 
associated  error  can  be  unnecessarily  large.  Instead  estimates  are  derived  from  ensemble  averages 
such  as  the  sample  covariance.  To  bridge  the  gap  between  these  two  approaches,  the  Fisher 
Information  for  the  sample  covariance  is  shown  to  be  equal  to  the  number  of  independent  and 
stationary  samples  times  the  Fisher  Information  for  a  single  sample.  Therefore,  there  are  no 
practical  limits  on  parameter  resolution  if  (1)  the  bound  for  a  single  sample  is  finite,  which  is 
generally  the  case  of  interest,  (2)  a  sufficiently  large  population  of  independent  samples  can  be 
found.  The  parameter  resolution  issue  then  becomes  one  of  determining  the  maximum  number  of 
such  samples.  This  number  is  set  by  physical  variables  that  do  not  appear  in  the  monochromatic  or 
instantaneous  measurement.  A  means  of  determining  this  number  from  the  temporal  coherence  of 
the  received  field  and  the  measurement  time  is  presented. 
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INTRODUCTION 

Recently  an  inconsistency  has  emerged  between  theo¬ 
retical  and  applied  ocean  acoustics.  A  number  of  parameter 
resolution  bounds  have  appeared  in  the  literature.  Many  of 
these  bounds  have  analytic  forms  that  do  not  allow  for  the 
possibility  of  ensemble  averaging,  including  some  bounds 
given  by  the  present  author,  and  therefore  predict  unrealisti¬ 
cally  poor  resolution  for  a  given  set  of  parameters. 1-5  This  is 
in  direct  contrast  to  current  practice  in  applied  ocean  acous¬ 
tics  where  ensemble  averaging  is  widely  used  to  lower  the 
variance  in  parameter  estimation. 

In  particular,  the  complex  sample  covariance  of  a  sensor 
array  is  generally  the  measurement  from  which  parameters 
are  estimated  in  applications  of  narrow-band  matched  field 
processing  and  beamforming.6  The  sample  covariance  is  ob¬ 
tained  by  ensemble  averaging  instantaneous  measurements 
of  a  sensor  array’s  spatial  covariance  for  a  fixed  narrow-band 
about  a  given  frequency.  While  the  averaging  is  usually  done 
by  sampling  in  the  time  domain,  the  result  is  equivalent  to 
averaging  the  same  number  of  independent  frequency  com¬ 
ponents  across  the  narrow-band  of  the  measurement.  As  an 
estimator  of  the  true  covariance,  the  sample  covariance  has 
obvious  advantages  over  an  instantaneous  measurement  be¬ 
cause  its  variance  can  be  linearly  reduced  with  the  number  of 
independent  and  stationary  samples  in  the  estimate. 

However,  many  theoretical  parameter  resolution  bounds 
for  underwater  acoustics  collapse  the  finite  bandwidth  of  the 
sample  covariance  to  a  single-frequency  or  monochromatic 
measure. 1-5  But  a  monochromatic  measurement  is  statisti¬ 
cally  similar  to  an  instantaneous  one  in  the  sense  that  each 
constitutes  only  a  single  independent  sample  of  a  temporally 
fluctuating  quantity.  This  is  because  the  time-bandwidth 
product  of  a  monochromatic  measurement,  like  that  of  an 
instantaneous  one,  is  by  definition  unity.  Clearly,  the  vari¬ 
ance  of  a  parameter  estimate  derived  from  a  single  sample 


may  be  unnecessarily  large  if  additional  independent  samples 
are  available.  But  it  is  implicit  in  ocean-acoustic  practice  that 
such  additional  samples  are  typically  available  because  the 
very  evidence  that  is  necessary  to  establish  the  validity  of  a 
statistical  approach  must  be  obtained  from  the  observation  of 
far  more  than  a  single  random  sample  of  the  measurement 
variable.  This  same  abundance  of  samples  can  then  be  used 
to  reduce  the  variance  in  a  parameter  estimate. 

To  bridge  the  gap  between  theory  and  practice,  the 
Fisher  information  for  a  measurement  of  the  complex  sample 
covariance  is  shown  to  be  equal  to  the  number  of  indepen¬ 
dent  and  stationary  samples  available  times  the  Fisher  infor¬ 
mation  for  a  single  sample.  Therefore,  there  are  no  practical 
limits  on  parameter  resolution  if  (1)  the  bound  for  a  single 
sample  is  finite,  which  is  generally  the  case  of  interest,  (2)  a 
sufficiently  large  population  of  such  samples  can  be  found. 

The  parameter  resolution  issue  then  becomes  one  of  de¬ 
termining  the  maximum  number  of  independent  and  identi¬ 
cally  distributed  samples  available.  However,  this  number  is 
set  by  physical  variables  that  do  not  appear  in  the  instanta¬ 
neous  or  monochromatic  measurement,  and  consequently,  is 
not  accounted  for  in  the  monochromatic  bounds  of  the 
present  ocean-acoustics  literature. 1-5  To  remedy  this  situa¬ 
tion,  a  means  of  determining  the  maximum  number  of  inde¬ 
pendent  samples  available  in  a  given  stationary  measurement 
period  is  presented. 

It  is  notable  that  in  their  matched  field  treatments  of 
broadband  signals.  Song7'8  as  well  as  Fawcett  and  Maranda9 
have  shown  that  the  optimal  position  resolution  of  a  target  in 
a  waveguide  can  generally  improve  by  increasing  the  dura¬ 
tion  of  the  received  signal,  if  other  parameters  such  as  the 
signal  bandwidth  and  noise  characteristics  are  held  fixed. 
This  is  consistent  with  well-established  radar  range- 
estimation  bounds  derived  for  targets  in  free  space.10  Appar¬ 
ently,  the  implication  of  these  results  has  not  been  fully  re- 
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alized  in  ocean-acoustic  modeling  of  narrow-band  signals. 
Therefore,  it  is  the  aim  of  the  present  paper  to  explicitly 
show  how  optimal  parameter  resolution  is  increased  by  av¬ 
eraging  independent  fluctuations  of  the  measurement  vari¬ 
able. 

In  Sec.  I,  a  general  expression  is  derived  for  the  Fisher 
information  matrix  contained  in  a  measurement  of  the  com¬ 
plex  sample  covariance.  This  is  done  without  specific  knowl¬ 
edge  of  the  probability  distribution  of  the  measurement.  Two 
practical  cases  are  examined  in  Secs.  II  and  III,  respectively. 
First,  the  entire  received  field  is  assumed  to  be  fully  random¬ 
ized.  Then,  the  signal  is  allowed  to  have  a  nonfluctuating 
component. 

In  the  first  case,  it  is  assumed  that  the  instantaneous 
fields  received  by  the  sensor  array  are  circular  complex 
Gaussian  random  (CCGR)  variables11  so  that  both  the  signal 
and  noise  components  fluctuate  over  temporal  samples.  This 
is  a  standard  assumption  implicit  in  the  derivation  of  many 
ocean  acoustics  parameter  bounds, 1-5  but  is  also  well 
founded  in  a  variety  of  other  fields  such  as  optics  and 
radar.11'12  It  is  additionally  assumed  that  measurements 
across  the  array  may  be  correlated  at  any  instant,  but  all 
discrete  samples  are  independent  and  stationary  over  time. 
The  Fisher  information  matrix  for  the  sample  covariance  is 
then  found  to  equal  the  product  of  the  Fisher  information 
matrix  for  an  instantaneous  measurement  of  the  field  across 
the  array  and  the  number  of  independent  and  stationary 
samples.  Consequently,  the  Fisher  information  matrix  of 
Refs.  1-5  is  found  to  only  be  valid  for  the  special  case  of  a 
single  monochromatic  measurement. 

Previously  derived  monochromatic  parameter  resolution 
bounds,1  5  therefore,  are  reduced  in  direct  proportion  to  the 
available  number  of  independent  and  stationary  samples 
when  finite  bandwidth  is  allowed.  An  expression  for  this 
number  is  given  in  terms  of  the  temporal  coherence  of  the 
received  field  and  the  measurement  time.  Appropriately  sub¬ 
stituting  this  expression  into  an  analytic  formula  recently  de¬ 
rived  by  Lee  et  alP  for  the  exact  distribution  of  a  discrete- 
sample-covariance  element  yields  a  new  approximate 
probability  distribution  for  an  element  of  the  continuous-time 
sample  covariance.  This  approximate  distribution  is  the  com¬ 
plex  counterpart  to  Rice’s  distribution  for  mean-square  noise 
current14  and  Mandel’s  distribution  for  the  intensity  fluctua¬ 
tions  of  a  photon  beam.15  In  optical  applications,  it  can  be 
used  to  describe  the  statistical  properties  of  finite-time- 
averaged  mutual  intensity11  under  the  assumption  of  cross- 
spectral  purity.16 

In  the  second  case,  it  is  assumed  that  the  field  received 
by  the  sensor  array  contains  a  deterministic  signal  compo¬ 
nent  amid  an  additive  CCGR  component.  It  is  shown  that  the 
Fisher  information  matrix  for  this  sample  covariance  is  dif¬ 
ferent  from  the  fully  random  case  but  is  directly  proportional 
to  the  number  of  independent  samples  available,  as  antici¬ 
pated.  Optimal  parameter  resolution,  therefore,  can  again  be 
significantly  improved  by  averaging  identically  distributed 
instantaneous  covariance  samples.  The  number  of  indepen¬ 
dent  samples  available  is  again  given  in  terms  of  the  tempo¬ 
ral  coherence  of  the  random  component  of  the  received  field 
and  the  measurement  time. 


Two  illustrative  examples  that  have  a  significant  impact 
on  the  interpretation  of  published  literature  are  given  in  Sec. 
IV.  These  are  for  matched  field  tomography  of  internal 
waves,  and  target  detection  using  ambient  noise  as  a  source 
of  opportunity. 

I.  FISHER  INFORMATION  FOR  THE  COMPLEX 
SAMPLE  COVARIANCE 

The  random  vector  *[n]  contains  the  instantaneous 
complex  fields  <J>j[n]  measured  by  sensors  i  =  l,2,3,...,/V<|,  at 
discrete  time  n.  To  permit  spatial  coherence  across  the  array 
at  any  instant  n,  it  is  assumed  that  any  sample  </>,[«]  may  be 
correlated  with  any  <f>j[n].  However,  all  field  samples  with 
differing  discrete-time  indexes  are  assumed  to  be  indepen¬ 
dent.  For  the  intended  application  in  interferometry,  all  fields 
are  assumed  to  occupy  a  narrow  band  about  the  same  carrier 
frequency  such  that  the  bandwidth  is  much  less  than  the 
propagation  speed  divided  by  the  array  aperture.  More  pre¬ 
cise  constraints  on  the  bandwidth  for  application  in  matched 
field  processing  are  discussed  in  Sec.  IV.  A  specific  exten¬ 
sion  to  broadband  applications  is  provided  for  a  fully  ran¬ 
domized  field  in  Sec.  II. 

The  complex  sample  covariance  is  then  defined  as 
1  N 

S[Af]=-2  <&[«]<!>>],  (1) 

Jy  n=  1 

where  S[A]  is  a  Hermitian  matrix  containing  the  elemental 
measurement  statistics  Stj[N],  for  ij =  1,2,3,..., N& ,  where 
H  denotes  complex  conjugate  transpose.  Here,  the  instanta¬ 
neous  mutual  intensity11  samples  <&[n]<&H[n]  are  assumed 
to  be  identically  distributed  over  time  index  n .  Practical  sce¬ 
narios  for  which  this  assumption  is  valid  are  investigated  in 
Secs.  II  and  III. 

The  complex  vector  W  of  dimension  N%  is  formed  by 
stacking  the  columns  of  S,  so  that  the  elements  of  W  are 
given  by  Wjj=Sjj[N],  for  i,j=l,2,3,...,N$,  where  the 
double  subscripted  vector  element  follows  notation  discussed 
in  the  Appendix  for  the  Kronecker  product.17  It  is  notewor¬ 
thy  that  in  this  context  Wr=W*.  The  statistics  of  W  are 
characterized  by  mean  vector  M=(W),  with  elements  M , 
and  covariance  matrix  C (N),  with  elements 

C,j,kl(N)  =  ( ( Wi}  ~MU)(  WkI  -  Mkl)  * ) , 

where  i,j,k,l=  1,2,3, .  The  dependence  of  the  covari¬ 
ance  C(N)  on  the  number  of  independent  and  stationary 
samples  N  is  explicitly  indicated  for  future  reference. 

The  vector  g  contains  the  parameters  gk  for 
k  =  1,2,3, . . .  ,Ng  .  These  parameters  are  assumed  to  be  nonsto¬ 
chastic  but  presumably  can  be  estimated  from  the  stochastic 
data  W. 

While  the  general  distribution  for  W  remains  unknown 
and  may  be  unrealizable  in  compact  form,18  it  is  assumed  to 
have  a  Gaussian  asymptote  for  large  N.  This  Gaussian  as¬ 
ymptote  can  be  written  as  a  conditional  probability,  for  W 
given  g,  by  the  relation 

exp{-i[W-M(g)]I’C-I(g;iV)[W-M(g)]} 

pwg- - 2 - ■ 

((2rr)lvi|C(g;V)|)1/2 

(2) 
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To  be  consistent  with  the  Kronecker  product  and  the  fact  that 
both  W  and  M  have  double  subscripts  and  contain  all  the 
elements  of  corresponding  Hermitian  matrices,  the  algebra  is 
akin  to  vector  ’  rather  than  complex  vector  ’  algebra. 
Specifically,  the  transpose  in  the  exponent  is  just  the  trans¬ 
pose  and  not  the  Hermitian  transpose,  and  the  factors  2  and 
\  ,  as  well  as  the  power  \  are  appropriately  designated  in  Eq. 
(2). 

For  estimation  of  the  parameters  g,  the  elements  of  the 
Fisher  information  matrix  J(g)  are  defined  by19 

d2 

Ju(g)  =  -E  — —  lnP(W|g)  .  (3) 

L dgidgj  \ 

Since  the  exact  form  of  the  distribution  for  R(W|g)  is  un¬ 
known,  it  seems  that  the  Fisher  information  matrix  for  the 
sample  covariance  cannot  be  determined  directly  from  Eq. 
(3).  However,  there  is  an  alternative  approach. 

It  is  well  established  that  if  g  is  a  function  of  the  Na 
parameters  in  vector  a  that  have  a  known  Fisher  information 
matrix  J(a),  the  Fisher  information  matrix  J(g)  can  be  deter¬ 
mined  from19’20 


X,  ,  ^  x-1,  X  d^T 

J(s(a))-^TJ  (a,^ST 


Suppose  the  parameter  a  is  the  mean  of  W  such  that  a=M. 
Then  because  P(Wja)  is  asymptotically  Gaussian  for  suffi¬ 
ciently  large  N,  it  is  also  well  established  that  the  inverse 
Fisher  information  matrix  J  1  ( M ) ,  in  an  estimate  of  the 
mean  M,  is  equal  to  the  covariance  matrix  C(A0  for  all 
A'3=  1 ,20  Therefore,  substituting  M  for  a  and  then  C(  N)  for 
J  1  ( M )  in  Eq.  (4)  yields  the  Fisher  information  matrix 
J(g(M))  for  all  1 . 

This  same  result  can  also  be  obtained  by  loosely  follow¬ 
ing  Fisher’s  use  of  the  central  limit  theorem  in  his  initial 
derivation21  of  what  later  became  known  as  Fisher  informa¬ 
tion.  First  note  that  C(A0=C(1)/A  must  hold  due  to  the 
independence  and  stationarity  of  the  N  samples  in  the  sample 
covariance  of  Eq.  (1).  Substituting  the  Gaussian  asymptote 
of  m\'|g)  into  the  right-hand  side  of  Eq.  (3)  then  yields  the 
well-known  result20 


C  ‘(1) 


*C(1) 

dgi 


c  ‘(1) 


<?C(1) 

dgj 


(9Mr 

+  N  —  C-‘(l) 
dgi 


<9M 

dgj’ 


(5) 

which  in  the  present  case  approaches  Jjj(g)  as  N  increases. 
Clearly,  only  the  term  involving  partial  derivatives  of  M 
with  respect  to  the  parameters  to  be  estimated  survives  as  N 
becomes  sufficiently  large.  But  the  resulting  asymptotic  ex¬ 
pression  for  J ijig)  is  valid  for  all  N'S=  1  because  the  Fisher 
information  for  the  sum  of  N  independent  and  identically 
distributed  joint  measurements  equals  N  times  the  Fisher  in¬ 
formation  for  a  single  such  measurement.22  For  example,  this 
is  why  J(M)=iVC_1(l)  must  hold. 

The  general  Fisher  information  matrix  J(g)  for  the 
sample  covariance  vector  W  can  then  be  written  element  by 
element  as 


c)Mt  ,  <?M 

Jij(g)=N  —  C- ‘(1)  — .  (6) 

1  dgi  dgj 

Again,  this  expression  is  exact  for  the  general  case  of  N 3*  1 . 
Further,  application  of  Eq.  (6)  is  not  limited  to  the  Hermitian 
sample  covariance  measurements  of  Eq.  (1),  but  is  valid  for 
any  ensemble  average  of  N  independent  and  identically  dis¬ 
tributed  samples  that  has  asymptotically  Gaussian  statistics 
for  N  sufficiently  large,  as  follows  from  the  referenced  works 
of  Fisher  and  Rao  for  real  measurement  vectors. 

The  direct  proportionality  to  sample  size  exhibited  in 
Eq.  (6)  may  seem  obvious  because  Fisher  information  is  ad- 
ditive  for  independent  samples  ’  ’  and  is  not  lost  when 
such  independent  identically  distributed  samples  are 
combined22  as  in  Eq.  (1).  However,  even  with  these  facts,  the 
preceding  analysis  side-stepped  the  difficulty  in  first  obtain¬ 
ing  an  expression  for  the  joint  distribution  of  the  instanta¬ 
neous  covariance  elements,  which  may  not  have  an  analyti¬ 
cally  realizable  form,  and  then  substituting  this  expression 
into  Eq.  (3). 

Equation  (6)  represents  the  Fisher  information  for  an 
arbitrary  population  of  N  independent  samples,  which  may 
be  collected  over  time  or  space,  so  long  as  the  statistics  are 
ergodic.11  If  the  statistics  are  not  ergodic  over  the  entire 
population,  the  expression  must  be  modified  so  as  to  add  an 
additional  term  of  the  form  of  Eq.  (6)  for  each  ergodic  sub¬ 
population  with  its  own  mean,  covariance  and  sample  size.  If 
the  measurements  are  restricted  to  temporal  samples  that  are 
not  stationary,  the  expression  must  again  be  modified  so  as  to 
add  an  additional  term  of  the  form  of  Eq.  (6)  for  each  sta¬ 
tionary  subpopulation.  The  merit  of  this  approach  becomes 
questionable  if  the  parameters  to  be  estimated  vary  signifi¬ 
cantly  during  what  was  presumed  to  be  a  stationary  measure¬ 
ment.  For  simplicity,  from  hereon  the  stochastic  component 
of  the  field  received  at  a  given  hydrophone  is  assumed  to  be 
stationary  in  the  sense  that  its  temporal  coherence  depends 
only  upon  the  time  difference  between  any  two  measure¬ 
ments. 

It  is  useful  to  make  the  dependence  on  N  explicit  in  the 
Fisher  information  matrix  J(g \N)  so  that  the  additive  prop¬ 
erty  of  the  information  in  independent  identically  distributed 
samples  is  clearly  exhibited  via  J(g;A)  =  AJ(g;l).  Then  the 
Cramer-Rao  lower  bound19  20  is  given  as 

£[(l,-«i)2]^[J“1(g;Af)]/i=^[J-1(g;D]//,  (7) 

which  unambiguously  shows  that  the  limiting  mean-square 
estimation  error  for  any  unbiased  estimate  gt  of  parameter  g , 
can  be  made  arbitrarily  small  if  (1)  a  sufficiently  large  popu¬ 
lation  N  of  independent  and  identically  distributed  samples 
can  be  found,  and  (2)  the  N=  1  single-sample  bound  is  finite. 

In  summary,  Eqs.  (6)  and  (7)  bound  the  minimum  error 
in  estimating  parameters  from  the  average  of  N  independent 
and  identically  distributed  measurement  samples,  and  exhibit 
a  well-known  inverse  dependence  on  sample  size  N. 19-22 
Some  specific  cases  are  now  considered. 
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II.  THE  RECEIVED  FIELD  IS  FULLY  RANDOMIZED 

A.  Fisher  information  for  the  complex  sample 
covariance 


It  is  assumed  that  the  instantaneous  field  samples  <f>j[n] 
are  CCGR  variables  regardless  of  the  relative  amplitude  be¬ 
tween  the  signal  and  noise  components.  It  is  again  assumed 
that  the  field  can  be  spatially  correlated  across  the  array  at 
any  instant,  and  that  all  field  measurements  with  differing 
discrete-time  indexes  are  independent. 

The  expectation  value  of  the  sample  covariance  has  el¬ 
ements  given  by 


1 

N 


N 


2 


(<M«]  <£*[«])> 


(8) 


where  M(,  is  independent  of  time.  By  use  of  CCGR  moment 
factoring,  1  the  covariance  of  the  sample  covariance  is  found 
to  have  the  elements 


Cqr,stW  =  E\ 


42  2  (</>?[n]0*[w])(^s[w]0*[m])* 

n=  1  m=  1 


—  MqrM* 
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The  Fisher  information  matrix  of  Eq.  (6)  can  be  written  as 

*[ML 

dgj 

(10) 

Following  the  Appendix,  the  identity  [CD(1)]9,.  st 
*]rf  is  employed.  By  also  noting  that 
[Mr]9,.  =  [M],.?  and  [M*_1]rf=[M~1]fr,  the  Fisher  informa¬ 
tion  matrix  can  be  written  as 
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Returning  again  to  more  familiar  double  subscripted  matrix 
notation,  let  the  expected  covariance  of  an  instantaneous 
sample  of  the  field  received  across  the  array  be  denoted  by 
the  matrix  K=(<1>[m]Oh[«]),  which  when  converted  to  a 
double  subscripted  vector  by  stacking  its  columns  yields  M. 
For  example,  the  relationship  between  the  elements  of  K  and 
M  is  given  by  Kij  =  M • .  Equation  (11)  can  then  be  written 
as 


f  ,  <?K  Dk1 

/,,(**)=*  ^K-'-K-'^J.  (.2, 

Because  the  signal  and  noise  are  assumed  to  be  fully  ran¬ 
domized,  absolute  phase  information  is  lost  between  inde¬ 
pendent  temporal  samples.  In  this  case,  there  is  no  difference 
between  the  Fisher  information  for  a  single  instantaneous 
sample  <!>[/;]  of  the  field  across  the  array  and  that  for  a 
single  instantaneous  sample  <I>[n]<l>H[n]  of  the  correspond¬ 
ing  mutual  intensities.  Therefore  Eq.  (12)  also  defines  the 


Fisher  information  matrix  for  N  independent  and  identically 
distributed  joint  measurements  of  the  instantaneous  field 
®[n]  across  the  array.  This  last  observation  provides  a  use¬ 
ful  confirmation  of  the  foregoing  derivation  because  it  en¬ 
ables  Eq.  (12)  to  be  immediately  obtained  from  the  well- 
known  form  of  the  Fisher  information  matrix  for  a  fully 
randomized  zero-mean  complex  Gaussian  measurement 
vector20  and  the  additive  property  of  Fisher  information  for 
independent  measurements.19  20,22 

The  Fisher  information  matrix  used  in  Refs.  1-5  for  a 
monochromatic  measurement  of  the  temporally  fluctuating 
field  received  by  a  hydrophone  array  is  only  valid  for  the 
special  case  when  N=  1  and  the  measurement  is  of  suffi¬ 
ciently  narrow  bandwidth.  As  noted  in  the  introduction,  how¬ 
ever,  a  single  sample  is  the  worst  possible  sample  popula¬ 
tion.  Therefore,  a  Cramer-Rao  lower  bound  derived  for 
N=  1  has  questionable  relevance  as  a  “fundamental”  bound 
on  parameter  resolution  because  it  gives  the  best  possible 
parameter  resolution  for  the  worst  possible  number  of 
samples.  The  ramifications  of  Eq.  (12)  on  specific  results  in 
the  literature  are  fleshed  out  in  Sec.  IV.  An  evaluation  of  the 
number  of  independent  samples  available  in  a  given  time 
period  under  the  CCGR  field  assumption  is  presented  next. 


B.  Sample  size  as  a  function  of  time  and  coherence 

An  expression  for  the  maximum  number  of  independent 
samples  available  in  a  stationary  measurement  period  is  now 
derived.  This  is  given  in  terms  of  the  temporal  coherence  of 
the  received  field  and  the  measurement  time.  The  general 
approach  is  analogous  to  that  used  by  Rice14  and  Mandel15 
for  related  problems.  In  loose  terms,  the  concept  is  to  deter¬ 
mine  the  number  of  times  the  received  field  is  expected  to 
fluctuate  independently  during  the  given  measurement  pe¬ 
riod.  This  is  achieved  by  inspection  of  the  signal-to-noise 
ratio  (SNR)  of  the  measurement.  Here  the  SNR  is  defined  as 
the  squared-mean  to  variance  ratio.  For  the  discretely 
sampled  case,  the  SNR  of  a  sample  covariance  element  is 


SNR{S,,[A]}  = 


I  Mi 


Cijjj(N) 


=N 


M: 


MuMjj 


(13) 


where  Mu  is  positive  semidefinite  and  equal  to  the  expected 
intensity  at  sensor  i.  Here  the  number  of  independent 
samples  N  is  equal  to  the  SNR  for  a  diagonal  element  of  the 
sample  covariance,  such  that  SNR{5,-([V]}  =  N.  This  is  be¬ 
cause  5,,[1]  has  an  expectation  value  that  equals  its  standard 
deviation  under  the  CCGR  field  assumption,  and  all  N 
samples  are  independent  and  identically  distributed. 

Analogously,  the  number  of  independent  samples  avail¬ 
able  in  a  continuous  time  measurement  of  Su  is  given  by  its 
SNR.  To  show  this,  the  sample  covariance  of  Eq.  (1)  can  be 
equivalently  written  as  a  continuous  temporal  average 


1  cm 

S  (T)=-  < (14) 
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where  the  received  fields  4>(f)  are  again  assumed  to  be  nar¬ 
row  band.  For  the  continuous  case,  the  SNR  of  •S',/7’)  is 
defined  as 


SNR{Sy(r)}  = 


</?y(r))2  +  </y(D)2 
2  2 
a~Ru(T)+a"iij(T) 


(15) 


where  Rij(T)=Re{Sij(T)},  Iu(T)=Im{Su(T)},  a\  <T)  is 

2 

the  variance  of  RjfT),  and  ,  i  >  is  the  variance  of  /;;(7).  It 
is  not  difficult  to  show  that  {Rij(T))=Re{Mij}  and 

</yCT))=Im{My}. 

Expressions  for  the  variances  can  also  be  obtained,  but 
with  more  difficulty.  First,  it  is  useful  to  employ  some  defi¬ 
nitions  from  statistical  optics.  The  complex  degree  of 
coherence11  is  defined  as 


Viji  t)  =  { t)  ip*  (t))/(M  hM  jj)  m,  (16) 

and  the  complex  coherence  factor1 1  is  defined  as 

Vij=7ij(0  )=Mij/(MiiMjj)m.  (17) 

Because  the  complex  degree  of  coherence  y/;(  t)  only  de¬ 
pends  on  the  time  difference  r,  the  fields  <b,(t)  must  be  sta¬ 
tionary  over  the  measurement  period  T  for  this  analysis  to  be 
valid. 

By  defining  the  normalized  cross-spectral  density 
3? ij(f )  as  the  Fourier  transform  of  the  complex  degree  of 
coherence  yy(r),  and  Qjj(f)  =  (MiiMjj)m3(?ij(f)  as  the  un¬ 
normalized  cross-spectral  density  which  is  the  Fourier  trans¬ 
form  of  the  mutual  coherence  function  ( <^>,(  t  +  t)  <f>*  (t)) ,  the 
expectation  value  of  the  sample  covariance  can  be  expressed 
as 


*0 

*  — 


Mu=  |  Qij(f)df. 


(18) 


This  spectral  representation  is  useful  for  computing  the 
Fisher  information  matrix  of  Eq.  (12)  with  state-of-the-art 
ocean-acoustic  propagation  models. 

For  illustrative  purposes,  it  is  now  assumed  that  the 
measurements  across  the  array  are  cross-spectrally  pure.16 
The  mathematical  expression  of  cross-spectral  purity  is 
yij(  t )  =  vt j  y(  r) ,  where  y(  r)  =  yi;(  r) ,  for  all  sensors  i  and  j . 
This  means  that  the  normalized  cross-spectral  density  mea¬ 
sured  between  any  two  spatially  separated  sensors  is  the 
same  as  the  normalized  spectral  density  measured  at  any  sen¬ 
sor  in  the  array.  While  this  condition  is  most  easily  satisfied 
for  narrow-band  measurements,16  which  are  most  commonly 
used  in  ocean-acoustic  interferometry,1'9  it  is  also  applicable 
to  broadband  processes. 

With  these  definitions,  expressions  for  the  variances  of 
the  real  and  imaginary  parts  of  the  sample  covariance  ele¬ 
ments  can  be  obtained  by  extending  some  results  presented 
by  Goodman, 1 1  in  his  discussion  of  the  properties  of  mutual 
intensity,  to  more  the  general  forms 
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where  R,j(T )  and  /(;(7’)  are  found  to  be  uncorrelated.  Here, 
/i  is  defined  by 

(20) 

and  the  triangle  function  is  defined  as 


A  ( r)  =  1  —  |  r| ,  for  |  r|  =£  1 , 

(21a) 

A(r)  =  0,  elsewhere. 

(21b) 

In  terms  of  the  spectral  density  -(S(f) ,  which  is  the  Fourier 
transform  of  yir),  a  useful  spectral  representation  for  /i  is 
given  by 


sin(  TvT{f-f')\ 

J 


df  df 


■  (22) 


By  appropriate  substitution  of  the  means  and  variances 
given  above,  Eq.  (15)  yields  the  desired  SNR  for  the  con¬ 
tinuous  measurement  case 

SNR{5(J(r)}  =  /r|  Vjj\2.  (23) 

Therefore,  the  number  of  independent  samples  available  in 
continuous  measurement  time  T  is  /i,  where  /i  need  not  be 
discrete  but  must  be  greater  than  or  equal  to  one,  as  is  evi¬ 
dent  by  inspection  of  Eqs.  (20)  and  (22).  It  is  important  to 
realize  that  if  the  assumption  of  cross-spectral  purity  cannot 
be  made,  the  number  of  independent  samples  SNR{5';,(7’)} 
would  not  necessarily  be  identical  across  all  sensors  i  of  the 
array,  but  would  be  given  by  Eqs.  (20)  or  (22)  with  yir)  and 
.y:(f)  replaced  by  y;,(r)  and  ,,(/ ),  respectively. 

The  continuous  sample  size  /i  can  also  be  interpreted  as 
the  time-bandwidth  product  of  the  cross-spectrally  pure  field 
received  by  the  array. 11,23  For  example,  by  defining  the  co¬ 
herence  time  scale  as 

r  oo  r  oo 

tc=  |y(r)|2c/r=  \5^\f)\2  df,  (24) 

J  —00  J  —00 

where  1  Itc  measures  the  bandwidth  of  field  fluctuations  over 
an  infinite  time  window,  either  Eq.  (20)  or  (22)  can  be  used 
to  show  that  as  T  becomes  much  greater  than  r,  ,  /i  ap¬ 
proaches  77 tc.  Specifically,  the  triangle  function  approaches 
a  constant  in  Eq.  (20)  and  the  sine-squared  function  ap¬ 
proaches  a  delta  function  in  Eq.  (22).  Therefore,  the  effect  of 
the  finite  measurement  window  on  the  bandwidth  of  the  re¬ 
ceived  field  B  =  /i/T  becomes  negligible  in  this  limit  so  that 
the  time-bandwidth  product  /i  grows  linearly  with  measure¬ 
ment  time  T. 

In  the  opposite  extreme,  inspection  of  Eq.  (20)  or  (22) 
indicates  that  as  tc  becomes  much  greater  than  T,  the  time- 
bandwidth  product  /i  approaches  unity  and  can  no  longer  be 
approximated  by  TItc.  This  is  because  the  coherence  time 
scale  of  Eq.  (24)  is  independent  of  the  temporal  window 
associated  with  a  specific  measurement.  Here,  the  complex 
degree  of  coherence  behaves  as  a  constant  in  Eq.  (20)  and  the 
spectral  density  behaves  as  a  delta  function  in  Eq.  (22).  This 
limiting  case  describes  both  the  monochromatic  and  the  in- 
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stantaneous  measurement,  for  which  the  number  of  indepen¬ 
dent  samples  must  be  unity,  as  is  well  known.11  With  regard 
to  the  introduction  of  this  paper,  the  foregoing  asymptotic 
analysis  provides  a  mathematical  basis  for  many  of  the 
points  raised  concerning  monochromatic  measures  and 
sample  size.  Specifically,  the  monochromatic  measurement 
is  the  subset  of  instantaneous  measurements  for  which  the 
spectral  density  -'/(f  )  not  only  behaves  as  a  delta  function  in 
Eq.  (22)  because  tc>T,  but  actually  is  a  delta  function 
.y(f)  =  S(f—f0)  in  the  sense  that  there  is  no  T  for  which 
tc9>T  is  not  true.  Therefore,  it  is  important  to  realize  that  a 
narrow-band  measurement  is  not  necessarily  monochro¬ 
matic,  but  a  monochromatic  model  is  often  sufficient  to  ap¬ 
proximate  the  interferometric  processes  involved  in  an  in¬ 
stantaneous  narrow-band  measurement.  However,  an 
instantaneous  broadband  measurement  must  generally  be  de¬ 
scribed  by  its  full  spectrum  ■':/ (f)  in  any  modeling  effort 
even  though  its  time-bandwidth  product  is  unity. 

It  is  useful  to  now  consider  the  influence  that  the  SNR  of 
Sjj(T)  has  on  estimators  found  in  ocean-acoustic  beamform¬ 
ing  and  matched  field  processing,  since  these  rely  upon  in¬ 
tersensor  coherence.  For  such  estimators,  Eq.  (23)  implies 
that  as  the  intersensor  measurements  become  less  correlated 
under  decreasing  vtj  | ,  longer  averaging  times  are  necessary 
to  achieve  the  same  resolution  for  a  given  parameter.  For 
practical  considerations,  the  temporal  augmentation  neces¬ 
sary  for  an  off-diagonal  element  at  i,j  to  attain  the  SNR  that 
a  diagonal  element  achieved  in  time  T  is  expected  to  be 
(Tl\vij\2-T),  when  T>tc. 

Finally,  the  minimum  sample  size  necessary,  but  not  al¬ 
ways  sufficient,  to  achieve  a  fixed  mean-square  error  toler¬ 
ance  in  a  parameter  estimate  can  be  equated  with  the  ratio  of 
the  Cramer-Rao  lower  bound  for  a  single  sample  to  the  error 
threshold.  This  minimum  value  then  sets  practical  require¬ 
ments  on  the  bandwidth  and  measurement  time  necessary  to 
attain  the  desired  resolution.  Illustrative  examples  along  this 
line  are  given  in  Sec.  IV. 


While  the  immediate  application  of  Fee  et  al. 1 3  was  for 
polarimetric  synthetic  aperture  radar  imaging  under  the 
CCGR  field  assumption,  their  exact  distributions  are  clearly 
useful  in  underwater  acoustic  interferometry.  Specifically,  in 
terms  of  its  complex  magnitude  aif  and  phase  f//(/ ,  the  dis¬ 
crete  sample  covariance  element  .S'(,[/V]  is  distributed  ac- 
cording  to 
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where  p,  is  the  magnitude  and  0U  is  the  phase  of  the  com- 
plex  coherence  factor  r>(/ ,  and  hij={M  UM  .  For  a  suffi¬ 
ciently  large  number  of  samples  N,  which  is  typically  greater 
than  ten  in  practice,  the  Gaussian  asymptote  for  is 

1  -T 

achieved.  The  magnitude  a,;  obeys  the  distribution 
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which  has  familiar  behavior  as  the  sensors  become  perfectly 
correlated,  where  p,.=  1,  and  uncorrelated,  where  p,  =0.  Fi- 
nally,  the  phase  (//,,  is  distributed  according  to 


C.  Approximate  distributions  for  finite-time-averaged 
mutual  intensity 

Under  the  CCGR  field  assumption,  the  approximate  dis¬ 
tribution  for  an  element  of  the  continuous-time  sample  cova¬ 
riance  Sjj(T)  can  be  obtained  by  combining  the  results  of  the 
last  section  with  an  exact  distribution  recently  derived  by 
Fee  et  al . 13  for  an  element  of  the  discrete  sample  covariance 
Sij[N],  This  approximate  distribution  is  the  complex  coun¬ 
terpart  to  Rice’s  distribution  for  mean-square  noise  current14 
and  Mandel’s  distribution  for  the  intensity  fluctuations  of  a 
photon  beam.15  In  acoustic,  optical,  and  radar  applications,  it 
can  be  used  to  describe  the  statistical  properties  of  finite- 
time-averaged  mutual  intensity11  under  the  assumption  of 
cross-spectral  purity.  The  idea  is  to  follow  the  analogous 
derivations  of  Rice  and  Mandel  by  substituting  the 
continuous-time  SNR  p  for  the  number  of  independent 
samples  N  in  the  exact  distribution  for  the  discrete  time  mea¬ 
surement  Sij[N]. 


T(jy+  1/2) ( 1  -  p2)NPij  cos(  <A-  0,j) 

^  2  1  -  (Pij  cos(  i/j-  0jj))2)N+ 1/2 

(l-p2.)^ 

+  0  11  F\N,  1;  l/2;(py  cos(  i/,-  %))2], 

Z  77 

(27) 

where  —  tt<  i/f,^  tt.  Here  standard  notation  is  used  for 
modified  Bessel  functions  7(),  K N .  , ,  gamma  function  T,  and 
Gauss’s  hypergeometric  function  F. 

Under  the  assumption  of  cross-spectral  purity  across  the 
array  aperture,  approximate  distributions  for  the  continuous¬ 
time  sample  covariance  elements  Sij(T)  of  Eq.  (14)  are  now 
obtained  by  substituting  the  time-bandwidth  product  p  of 
Eqs.  (20)  and  (22)  for  N  in  Eqs.  (25)-(27).  The  resulting 
distributions  approximately  describe  the  statistics  of  finite- 
time-averaged  mutual  intensity,  and  so  extend  Mandel’ s  ap- 
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proximate  distributions  for  the  diagonal  terms  Sti(T) 1 123  de¬ 
veloped  in  his  analysis  of  photon  beam  fluctuations.15  How¬ 
ever,  barring  the  trivial  case  in  which  the  elements  are  un¬ 
correlated,  the  joint  distribution  for  all  elements  of  the 
sample  covariance  apparently  remains  unrealized  in  analytic 
form,  as  noted  in  Sec.  I,  regardless  of  whether  the  sampling 
is  continuous  or  discrete. 

D.  Relationship  to  broadband  performance  bounds 

Besides  Refs.  7-9  there  are  a  number  of  papers  in  the 
literature  that  examine  ocean-acoustic  parameter  resolution 
bounds  for  broadband  measurements  on  a  hydrophone  array. 
For  example.  Refs.  24-26  provide  expressions  for  the  corre¬ 
sponding  Fisher  information  matrices  in  terms  of  the  wave¬ 
guide  Green  function.  However,  the  spectral  summations  and 
integrals  used  in  these  expressions  do  not  include  the  convo¬ 
luted  effects  of  finite  time  windows  nor  do  they  include  any 
description  of  temporal  sampling.  Exclusion  of  the  former 
effect  limits  the  validity  of  these  expressions  to  the  asymp¬ 
totic  regime  where  the  measurement’s  temporal  window  is 
so  much  greater  than  the  coherence  time  scale  that  the  asso¬ 
ciated  sample  size  must  be  statistically  large.  While  no  cor¬ 
responding  sample  size  is  explicitly  given  to  scale  the  Fisher 
information  matrix  in  these  references,  this  factor  may  be 
implicitly  defined.  Such  a  scaling  is  necessary  according  to 
Whittle’s  asymptotic  expression  for  the  Fisher  information 
matrix  for  a  set  of  time  series  measurements  with  large 
sample  size.27  However,  a  danger  in  such  asymptotic  ap¬ 
proximations  is  that  their  validity  may  inherently  be  limited 
to  sample  populations  of  an  enormity  that  cannot  be  achieved 
in  practice. 

While  the  Fisher  information  matrices  derived  in  Refs. 
28  and  29  also  neglect  the  effects  of  finite  time  windowing, 
they  do  so  by  direct  use  of  Whittle’s  approximate  Fisher 
information  matrix  and  explicitly  include  the  number  of  in¬ 
dependent  samples  available  in  the  measurement,  under  the 
assumption  that  this  number  is  large.  However,  the  acoustic 
ray-multipath  model  of  Refs.  28  and  29  cannot  be  generally 
reconciled  with  the  full-field  propagation  models  that  are 
presently  necessary  for  realistic  ocean-acoustic 
simulations. 1-9  Therefore  the  Fisher  information  matrices  de¬ 
rived  in  Refs.  28  and  29  are  of  limited  practical  use  in  ocean 
acoustics. 

The  present  analysis  stands  apart  from  these  studies  in 
that  (1)  it  is  not  confined  to  the  asymptotic  regime  but  is 
valid  for  the  general  case  of  N'S?  1 ,  and  (2)  it  can  be  easily 
implemented  with  the  requisite  full-field  propagation  models 
that  are  currently  state-of-the-art  in  ocean  acoustics.  While 
the  focus  is  on  narrow-band  interferometry,  all  of  the  expres¬ 
sions  derived  for  the  fully  randomized  field  case  of  Sec.  II 
are  directly  applicable  to  broadband  problems.  For  example, 
Eq.  (12)  is  not  limited  to  narrow-band  processes  but  also 
defines  the  Fisher  information  matrix  for  joint  measurements 
of  a  stationary  broadband  CCGR  field,  where  K  is  the  ex¬ 
pected  covariance  of  the  instantaneous  field  across  the  array, 
N  is  the  number  of  times  an  instantaneous  field  across  the 
array  is  measured  independently,  Eqs.  (20)  and  (22)  deter¬ 
mine  the  number  of  such  measurements  available  in  a  given 
period,  and  Eq.  (18)  gives  a  spectral  representation  for  the 


elements  of  K,  where  Kjj  =  M  • .  Finally,  it  is  noteworthy  that 
for  sufficiently  large  N,  the  general  expression  derived  here 
in  Eq.  (12)  is  consistent  with  Whittle’s  asymptotic  formula. 


III.  THE  RECEIVED  FIELD  HAS  A  DETERMINISTIC 
COMPONENT 

A.  Cramer-Rao  bounds  for  the  complex  sample 
covariance 


Let  the  instantaneous  field  measurements  <f>j[n]  contain 
a  deterministic  component  <fiD  [n]  superposed  with  a  CCGR 
component  <f>G  [«],  such  that  </>,[«]  =  4>D .[n]  +  cf)c  [n], 
which  is  equivalently  described  by  the  vector  equation 
<1>[h]  =<I>D[n] +<I>g[h] .  The  random  part  of  the  field 
4>c  [h]  may  contain  a  stochastic  signal  as  well  as  noise.  The 
deterministic  part  of  the  field  <f>D  [n]  is  assumed  to  occupy  a 
sufficiently  narrow-band  about  the  carrier  frequency  that 
4>d .[«]  =  4>d4>d is  time-invariant  during  collection 
of  the  N  samples.  It  is  again  assumed  that  the  field  can  be 
spatially  correlated  across  the  array  at  any  instant,  and  that 
all  field  measurements  with  differing  discrete-time  indexes 
are  independent. 

The  expectation  value  of  the  sample  covariance  has  el¬ 
ements  given  by 


1 

N 


N 


2 

n=  1 


=  (f>D  .(/>%.  +  ((f>Cj  <ko)  =  4>ol(f>Dj  +  Mcif  (28) 

where  each  element  contains  a  cross  term  4>d4>%.  for  the 
deterministic  field,  and  another  time-invariant  term  Mn  for 

ij 

the  random  field.  By  use  of  CCGR  moment  factoring,11  it 
can  be  shown  that  the  covariance  of  the  sample  covariance 
has  elements  given  by 


Cqr,sAN)=E 
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(09[n](/>*[n])(</>j[m]</>f[w])* 


—  M  M* 

1V1  qrlvl  st 

4  (MoJ*nr4>D+  tDqKM*G+MGqM*r),  (29a) 


which  becomes 

1 

Cqr,st(N)  =  -  (MGqs<l>*r(f>D+MqsM*ri),  (29b) 

by  application  of  Eq.  (28),  so  that  Cqr  st{N )  =  Cqr  st(  1  )/N  as 
anticipated. 

Therefore  Eq.  (10)  yields  the  desired  Fisher  information 
matrix  upon  substitution  of  the  mean  defined  in  Eq.  (28),  and 
the  covariance  defined  in  Eq.  (29a)  or  (29b),  where 
[C~1(A0]?l%,sf=iV[C  -1(l)]9r ,st  •  Clearly  this  Fisher  informa¬ 
tion  matrix  is  directly  proportional  to  ;V,  and  therefore  the 
Cramer-Rao  lower  bound  for  any  parameter  with  a  finite 
bound  for  N=  1  must  decrease  as  N  increases. 
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A  more  compact  expression  for  the  Cramer-Rao  lower 
bound 


c)Mt  ,  dM^' 
- C“‘(l) - 

dg  dg 


(30) 


is  easily  obtained  from  Eqs.  (6)  and  (7).  With  the  definitions 

[S ],•/  =  4>o  (f>f)  ,  and  =  Mc  ,  the  covariance  of  Eq. 

J  i  j  J  U 

(29)  may  be  written  as 


C  (AO 


1 

N 


(MGHr+  MMq). 


(31) 


This  leads  to  the  explicit  bound 


£[(i,-^/)2] 


{MGHr+MMG}-1 


d  M 

dg. 


which  may  be  computationally  expensive  to  implement  be¬ 
cause  it  requires  two  matrix  inversions.  It  is  noteworthy  that 
while  the  matrix  oM Idg  is  not  invertible  in  general,  the 
bound  may  always  be  written  in  the  form 


E[(g-gi)2~i^^ 


dg  dg1 

—  C(  1 )  — 
dM  v  ;  <?M 


(33) 


if  the  matrix  dg/dM  exists30  as  it  must  when  the  problem  is 
properly  constrained.  More  specifically,  this  bound  is  given 
by 


E[(g-gi)2]^^ 


{MGHr+MMG} 
dM  x  c  Gf  dM 


(34) 


This  form  requires  no  matrix  inversion,  but  it  may  be  less 
plausible  to  implement  than  Eq.  (32)  because  the  matrix 
dg/dM  is  usually  more  difficult  to  determine  than  dM! dg  for 
applications  in  ocean-acoustic  interferometry. 

Due  to  its  temporal  invariance  over  the  collection  of  the 
N  measurement  samples,  the  deterministic  component  H  is 
effectively  monochromatic.  By  defining  Qa.if)  as  the  Fou¬ 
rier  transform  of  the  mutual  coherence  function 


{4>c{t  +  T)4>%{t)), 

where  cf>G  (t)  is  a  continuous-time  measurement  of  the  ran¬ 
dom  field  component,  the  expected  covariance  MG  of  the 
random  field  component  can  be  expressed  spectrally  as 

r  co 

Mc  =  QG(f)df,  (35) 

J  J  —  00  J 


in  accord  with  the  analysis  of  Sec.  II  B.  This  representation 
of  Mc  is  convenient  for  numerical  implementation  of  the 
bound  given  in  Eqs.  (32)  and  (34)  because  many  state-of-the- 
art  ocean-acoustic  propagation  models  operate  in  the  spectral 
domain. 

It  is  noteworthy  that  when  the  deterministic  component 
vanishes,  so  that  H=0,  M  and  Mc  are  identical.  The  bound 
can  then  also  be  obtained  from  the  Fisher  information  matrix 
of  Eq.  (12),  for  the  fully  randomized  field  case,  and  is  there¬ 
fore  directly  applicable  to  broadband  measurements.  Con¬ 
versely,  if  the  random  component  of  the  field  vanishes  so 
that  Mg  =  0,  the  covariance  C (N)  is  zero.  In  this  case,  the 


bound  is  also  zero,  corresponding  to  perfect  resolution  for 
any  parameters  that  can  be  fully  described  by  their  depen¬ 
dence  on  the  mean  vector  M  and  known  constants. 

B.  Sample  size  as  a  function  of  time  and  coherence 

The  number  of  independent  samples  available  in  a  sta¬ 
tionary  measurement  period  is  given  by  the  time-bandwidth 
product  of  the  received  field’s  random  component.  This 
number  can  be  obtained  by  substituting 

((f>Cj{t+r)  <fiG{t))/M c.. 

for  yir)  in  Eq.  (20),  under  the  assumption  that  the  random 
component  is  cross-spectrally  pure  over  the  array  aperture. 
As  of  yet,  mutual-intensity  distributions  analogous  to  those 
of  Sec.  II  C  have  apparently  not  been  derived  for  the  case  of 
a  deterministic  signal  in  a  CCGR  field. 

It  is  significant  that  the  differing  assumptions  about  the 
received  field  made  in  Secs.  II  and  III  have  lead  to  differing 
Fisher  information  matrices.  Clearly,  the  statistics  of  the 
problem  at  hand  must  be  evaluated  before  bounds  can  be 
assigned.  Otherwise,  such  arbitrary  considerations  as  the 
mathematical  convenience  of  a  particular  formulation  of  the 
Fisher  information  begin  to  motivate  the  theoretical  analysis. 
It  is  equally  significant  that  these  results  are  unified  by  the 
fact  that  Fisher  information  increases  linearly  with  the  num¬ 
ber  of  independent  samples  used  to  construct  the  sample  co- 
variance,  which  is  a  natural  consequence  of  the  more  general 
formulation  given  in  Sec.  I  that  follows  standard  statistical 
methods. 

IV.  ILLUSTRATIVE  EXAMPLES 

The  Fisher  information  matrix  derived  in  Sec.  II  can  be 
directly  applied  to  appropriately  rescale  monochromatic 
bounds  that  already  exist  in  the  literature.  For  example,  the 
resolution  bounds  given  in  Refs.  1-5  can  be  made  more 
practical  by  dividing  the  limiting  minimum  mean-square  es¬ 
timation  error  presented  by  the  appropriate  number  of  inde¬ 
pendent  and  stationary  samples.  This  number  can  be  deter¬ 
mined  by  Eq.  (20)  or  (22)  if  the  temporal  coherence  of  the 
received  field  and  the  stationary  time  of  the  physical  pro¬ 
cesses  causing  the  fluctuations  are  known. 

Knowledge  of  the  maximum  allowable  sample  size  is 
most  critical  when  monochromatic  bounds  suggest  inad¬ 
equate  resolution  for  a  given  parameter.  Two  narrow-band 
examples  from  the  literature  are  presented.  In  both  examples, 
the  entire  received  field  is  assumed  to  obey  CCGR  statistics, 
and  the  original  Fisher  information  matrix  from  the  literature 
corresponds  to  that  given  by  Eq.  (12)  for  the  single  instanta¬ 
neous  measurement  case  of  N=  1,  where  the  measurement 
band  about  central  frequency  /0  is  sufficiently  narrow  to  ap¬ 
proximate  the  expected  instantaneous  covariance  K  of  Eq. 
(18),  where  Kjj=Mjj,  with  a  monochromatic  model  spec¬ 
trum  for  which  .^(f)  =  S(f— f0). 

A.  Matched-field  tomography  of  internal  waves 

First  consider  Ref.  5,  where  internal-wave  parameters 
are  to  be  estimated  by  full-field  acoustic  tomography  with  a 
controlled  source  signal  of  finite  duration  and  bandwidth.  To 
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resolve  a  long  internal  wave  to  roughly  an  acoustic  wave¬ 
length  \acou,  a  measurement  time  T  much  less  than  an  inter¬ 
nal  wave  period  rinl  is  generally  required,  under  the  assump¬ 
tion  that  \acou  is  much  less  than  the  internal  wavelength. 
Specifically,  the  maximum  measurement  time  rmax  should  be 
no  greater  than  Xacou/cnu.  For  a  typical  two-layer  shallow 
water  waveguide  with  density  difference  to  density  ratio 
Ap/p=10  3,  upper  layer  thickness  H=  20  m,  and  gravita¬ 
tional  acceleration  g,  the  long  internal  wave  speed 

Cint=(t?^p/p)1/2  (36) 

is  roughly  0.4  m/s.  To  obtain  more  than  a  single  sample,  the 
minimum  bandwidth  should  be  greater  than  cin/Xacou.  The 
restriction  that  the  processing  be  narrow  band  then  sets  the 
maximum  bandwidth  // max  and  sample  size  Nmax=7’max.Bmax . 
In  narrow-band  beamforming,  B  max  must  be  much  less  than 
the  propagation  speed  divided  by  the  array  aperture.  How¬ 
ever,  in  narrow-band  matched-field  processing,  the  propaga¬ 
tion  range  must  also  be  considered  to  insure  the  acoustic 
interference  structure  on  the  array  does  not  vary  significantly 
over  the  measurement  band.  In  this  case,  the  maximum  band¬ 
width  must  typically  be  determined  by  numerical  modeling. 

For  the  500-Hz  carrier  frequency  of  Ref.  5,  B max  must  be 
greater  than  1/Trnax=cint/A.acou=0.13  Hz  to  obtain  more  than 
a  single  sample.  But  parabolic  equation  modeling  shows  that 
B  max  can  be  much  greater  than  0.13  Hz  and  still  be  sufficient 
for  narrow-band  matched-field  inversion  over  the  ranges  and 
waveguide  discussed  in  Ref.  5,  so  that  A'nlax  can  greatly  ex¬ 
ceed  unity.  Therefore,  the  monochromatic  parameter  resolu¬ 
tion  bounds  presented  in  Ref.  5  may  be  significantly  im¬ 
proved  by  time  or  frequency  averaging.  Further,  the 
observation  of  that  reference  that  certain  parameters  may  not 
be  adequately  resolved,  even  if  the  number  of  spatial 
samples  with  unique  expectation  values  exceeds  the  number 
of  parameters  to  be  estimated,  is  not  necessarily  true  for  the 
specific  examples  given  therein  if  A'max  is  permitted  to  be 
greater  than  unity  in  accord  with  the  foregoing  analysis. 

Finally,  it  is  noteworthy  that  the  CCGR  fields  assumed 
for  this  example  may  in  reality  be  caused  by  such  diverse 
mechanisms  as  high-frequency  gravity-wave  induced  fluc¬ 
tuation  of  the  waveguide’s  surface  boundary  or  refractive 
index,  motion  of  the  source  or  receiver,  or  an  incoherent 
source. 

B.  Target  detection  using  ambient  noise  as  a  source 
of  opportunity 

Now  consider  Ref.  3,  where  an  object  submerged  in  a 
waveguide  is  to  be  detected  using  ambient  noise  as  a  source 
of  opportunity.  In  that  reference,  monochromatic  bounds  are 
derived  for  the  range,  cross-range  and  depth  resolution  of  the 
object.  But  the  limiting  resolutions  given  for  cross-range  and 
depth  are  extremely  poor  even  at  very  short  ranges  from  the 
object.  In  the  present  context,  those  mean-square  error 
bounds  become  more  realistic  when  they  are  divided  by  the 
maximum  allowable  sample  size  Amax .  For  example,  it  may 
be  possible  to  resolve  the  object  if  ;Vmax  is  sufficiently  large 
that  it  equals  the  corresponding  monochromatic  bound  of 
Ref.  3  divided  by  the  square  of  the  object  scale.  However, 
even  within  short  ranges  of  a  few  object  lengths,  this  ratio 


requires  /Vmax  to  be  greater  than  104.  When  Amax  is  equated 
with  the  time-bandwidth  product  of  the  received  signal,  it 
becomes  apparent  that  no  combination  of  temporal  coher¬ 
ence  and  stationarity  for  the  CCGR  ocean  surface  noise  pro¬ 
cess  will  make  that  detection  scheme  practical  beyond  a  few 
object  lengths.  The  conclusion  may  not  be  the  same  for 
highly  directional  noise,  however,  as  is  noted  in  Ref.  31. 

Therefore,  it  is  not  only  the  Fisher  information  for  a 
monochromatic  or  instantaneous  measurement,  but  also  the 
coherence  and  stationary  time  scales  of  the  random  compo¬ 
nent  of  the  received  field  that  are  essential  in  determining 
parameter  resolution  bounds  in  these  examples.  It  is  antici¬ 
pated  that  this  will  be  the  case  in  many  other  situations  of 
practical  interest  in  ocean  acoustics,  as  no  doubt  many  re¬ 
searchers  are  already  aware. 

V.  DISCUSSION  AND  CONCLUSIONS 

It  is  currently  common  practice  in  theoretical  ocean 
acoustics  to  derive  parameter  resolution  bounds  that  are  ex¬ 
clusively  for  a  monochromatic  measurement  of  the  fluctuat¬ 
ing  field  received  by  a  hydrophone  array.  Often  the  intent  of 
these  derivations  is  to  compare  the  resulting  bounds  with  the 
resolution  of  numerical  simulations  obtained  from  single¬ 
frequency  models  on  a  digital  computer.  It  is  only  natural 
that  such  comparisons  should  turn  out  favorably.  This  is  be¬ 
cause  the  single -frequency  model  inherently  assumes  zero 
bandwidth  which  is  equivalent  to  the  absence  of  temporal 
uncertainty — a  perpetually  instantaneous  field  structure.  A 
physically  relevant  analogy  to  this  situation  is  the  static  but 
random  interference  pattern  or  speckle  found  in  the  reflection 
of  a  coherent  laser  beam  from  a  smooth  surface  with 
wavelength-scale  roughness.32  In  both  of  these  cases,  the 
randomness  is  only  exhibited  by  a  spatial  sampling.  And  in 
both  cases  the  resulting  intensity  distribution  is  exponential 
with  a  standard  deviation  equal  to  the  mean,  implying  a  stan¬ 
dard  deviation  for  the  intensity  level  equal  to  Dyer’s  5.6 
dB,33  if  (1)  the  fields  are  sufficiently  randomized  to  be  sta¬ 
tionary  CCGR  variables,  and  (2)  the  spatial  aperture  is  less 
than  or  equal  to  the  speckle  scale. 

However,  once  finite  bandwidth  is  introduced  by  either 
(1)  an  incoherent  source,  (2)  relative  motion  between  source 
and  receiver,  (3)  relative  motion  between  scatterer  and 
source  or  receiver,  (4)  fluctuation  of  the  scatterer  or  wave¬ 
guide  boundary,  (5)  medium  scintillation,  or  (6)  the  simple 
introduction  of  noise,  the  randomness  takes  on  a  temporal 
dependence.  Averaging  the  intensity  in  time  at  any  point  in 
space  reduces  the  variance  by  the  number  of  independent  and 
stationary  fluctuations  N.  Similarly,  the  intensity-level  stan¬ 
dard  deviation  for  CCGR  fields  is  approximated  by  the  de¬ 
caying  quantity  (10  log  e)/Nul  dB  with  increasing  accuracy 
as  N  increases.23’34  In  this  context,  a  monochromatic  mea¬ 
surement  is  like  an  instantaneous  one  in  that  it  comprises  a 
single  sample  of  a  temporally  fluctuating  quantity.  But  a 
single  sample  constitutes  the  least  favorable  sample  popula¬ 
tion  as  far  as  parameter  estimation  is  concerned  because  it 
yields  the  largest  variance.  This  is  well  known  in  applied 
ocean  acoustics  where  parameter  estimates  are  typically  de¬ 
rived  from  ensemble  averages  such  as  the  complex  sample 
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covariance  which  is  widely  used  in  narrow-band  matched 
field  processing  and  beamforming. 

To  bridge  the  gap  between  these  current  theoretical  and 
applied  approaches,  the  Fisher  information  for  the  sample 
covariance  is  shown  to  be  equal  to  the  number  of  indepen¬ 
dent  and  stationary  samples  available  times  the  Fisher  infor¬ 
mation  for  a  single  sample.  Therefore,  there  are  no  practical 
limits  on  parameter  resolution  if  (1)  the  bound  for  a  single 
sample  is  finite,  which  is  generally  the  case  of  interest,  (2)  a 
sufficiently  large  population  of  independent  samples  can  be 
found.  The  parameter  resolution  issue  then  becomes  one  of 
determining  the  maximum  number  of  such  samples. 

A  means  of  determining  the  number  of  independent 
samples  available  in  a  given  measurement  period  is  pre¬ 
sented.  Some  specific  analytical  results  are  given  in  terms  of 
the  temporal  coherence  of  the  received  field,  which  simplify 
greatly  under  the  assumption  of  cross-spectral  purity  in  spa¬ 
tial  measurements  across  the  hydrophone  array.  This  ap¬ 
proach  has  also  led  to  a  means  of  applying  recently  derived 
probability  distributions  for  the  elements  of  the  discrete 
sample  covariance  to  approximate  similar  distributions  for 
finite-time-averaged  mutual  intensity. 

As  a  final  point  of  discussion,  the  Cramer-Rao  lower 
bound  is  typically  used  to  evaluate  the  performance  of  pa¬ 
rameter  estimators.  From  the  perspective  of  classical  estima¬ 
tion  theory, 19-22  optimal  estimators  are  unbiased  and  attain 
this  lower  bound.  While  optimal  estimators  are  only  certain 
to  be  found  if  the  parameters  to  be  estimated  are  linearly 
related  to  the  expectation  value  of  the  measurements,  maxi¬ 
mum  likelihood  estimators  are  asymptotically  optimal  for  a 
sufficiently  large  number  of  independent  measurement 
samples  even  if  the  relationship  between  the  parameter  set 
and  measurement  expectation  value  is  nonlinear. 19-22  This 
last  fact  is  of  particular  significance  in  ocean-acoustic  inter¬ 
ferometry  because  the  nonlinear  matched  field  estimators 
generally  required  may  only  attain  optimality  in  the  asymp¬ 
totic  limit  of  a  large  sample  population.  This  highlights  the 
significance  of  not  only  including  the  sample  size  as  a  vari¬ 
able  but  also  of  determining  its  maximum  allowable  value 
when  establishing  bounds  on  ocean-acoustic  parameter  reso¬ 
lution. 
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APPENDIX:  KRONECKER  PRODUCTS  AND  THEIR 
INVERSES 

From  Ref.  17,  the  Kronecker  product  of  the  double  sub¬ 
scripted  vectors  A  and  B  is  given  by 

[A®B],riJt=[A],J[B]rf.  (Al) 

When  the  inverses  A-1  and  B-1  exist,  then 

(A®B)”1  =  A-1«>B-1,  (A2) 

so  that 

[(A®B)-1],r,rt  =  [A-I]„[B-1]rt.  (A3) 


Therefore,  Eq.  (9)  indicates  that 
1 

C(A0=-M<g>M*.  (A4) 

Application  of  Eq.  (A3)  then  yields 

[C“1(iV)],r,,t=^[M-1],J[M*-1]rf,  (A5) 

under  the  assumptions  of  Sec.  II  A. 

For  notational  convenience,  the  definition  AB=A®B  is 
used  for  the  Kronecker  product  in  the  main  text. 
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